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In this letter we study finite temperature properties of spin-1/2 interacting bosons with spin- 
orbit coupling in two dimensions. When the ground state has stripe order, we show that thermal 
fluctuations will first melt the stripe order and lead to a superfiuid of boson pairs if the spin-orbit 
coupling is isotropic or nearly isotropic. Such a phase supports fractionalized quantum vortices. 
The Kosterlize-Thouless transition from superfiuid to normal state is driven by proliferation of half 
vortices. When the ground state is a plane wave state, the transition to normal state is driven by 
conventional Kosterlize-Thouless transition. However, the critical temperature will drop to zero for 
isotropic spin-orbit coupling. 



Spin-orbit coupling (SOC) plays an important role in a 
wide range of quantum systems from condensed matter, 
atomic physics to nuclear physics. For example, recent 
studies of electronic systems with strong SOC have led to 
the discovery of topological insulators [T] with potential 
applications to new quantum devices. Recently, bosons 
with SOC have also begun to attract considerable atten- 
tions [2H6], for such a coupling can now be generated in 
cold atom systems by engineering atom-light interaction 
[7]. Very recently, a pioneer experiment in NIST has al- 
ready achieved boson condensate with a restricted class of 
SOC E). Mean-field studies of the ground state have re- 
vealed two different phases, i.e. the plane wave phase and 
the stripe phase [SHSl |8] . However, the fluctuation effects 
on top of mean-field description and the finite tempera- 
ture properties have not been studied yet. In this letter 
we show that finite temperature thermal fluctuations will 
indeed lead to new phase and phase transitions. 

In this work we consider spin-1/2 interacting bosons 
with SOC in two dimensions, and present a global phase 
diagram of this system in terms of the interaction param- 
eters, temperature and the degree of anisotropy of SOC. 
The main results are summarized as follows: 

(1) For stripe phase with isotropic or nearly isotropic 
SOC, the system undergoes two phase transitions as tem- 
peratures increases. Stripe order melts at the lower tran- 
sition temperature and superfluidity disappears at the 
higher one. Between these two transitions, bosons form 
pairs and these boson pairs exhibit superfluidity. In this 
unusual boson paired phase, vortices are fractionalized. 

(2) For the stripe phase with anisotropic SOC, the sys- 
tem undergoes a direct Kosterlize-Thouless (KT) transi- 
tion to normal state. However, since the energy of a half 
vortex bound to a single dislocation is lower than that 
of a usual vortex, the KT transition is always driven by 
proliferation of half vortices. 

(3) For the plane wave phase, the system undergoes 
a conventional KT transition to normal phase driven by 
proliferation of usual vortices. However, the KT temper- 
ature drops to zero for isotropic SOC. 

We note such a boson paired phase and unconventional 
KT transition have also been proposed in other systems 



with different microscopic physics, such as stripe phase 
in high-Tc superconductor materials |9|, "Fulde-Ferrell- 
Larkin-Ovchinnikov" state [10] of superconductor in Zee- 
man field and spin-1 bosons [llj . however, it has not been 
realized experimentally yet. Ultracold bosons with SOC 
provide a new route toward this interesting physics, and 
may shed lights on understanding other systems. 
The single particle Hamiltonian Ho with SOC is 
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are Pauli matrices. Without loss of generality, we can 
define — k and with ^ 77 ^ 1. In this case, 
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with their wave function given by (p± 
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For interactions, generally speaking it can be written as 
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In Ref. W we considered the mean-field ground state of 
a simplified situation with goo = gii = go-: where "Hint 
can be rewritten as 
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where iV = np + ni, Sz = tlq — hi, and g = goi — ffoo- 
We found two distinct phases for g > or < 0. For 
5 < 0, all bosons condense into a spontaneously cho- 
sen single plane wave state, either ip^ or (/?_. This state 
breaks time-reversal symmetry. It is named "plane wave 
phase". For 5 > 0, bosons condense into a superposition 
state of these two plane wave states. Consequently, this 
state exhibits stripe order of spin density, and is named 
"stripe phase". Later studies show that this conclusion 
will not change qualitatively even if one considers more 
complicated interactions |5]. These two phases are anal- 
ogous to "Fulde-Ferrell" or "Larkin-Ovchinnikov" states 
proposed for superconductor in a Zeeman field. 
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FIG. 1: The amplitude of the order parameter with (a) no 
topological defect, (b) a single vortex, (c) double dislocations 
and (d) a half vortex bound to a dislocation. The red arrow 
indicates that the superfluid phase changes 2n around a single 
vortex (b) and n around a half vortex bound with a dislocation 
(d). 



For a typical situation in cold atoms, go S> \g\. Thus, 
we take the approximation that only spin and phase fluc- 
tuations are considered in the low-energy theory and the 
density fluctuations are ignored. First, let us consider 
the stripe phase, and we introduce the ansatz as 
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9 = {9+ + 9^)/2 and u = {9+ - 9^)/2 describe two 
low-lying degrees of freedom. 9 is the total phase which 
corresponds to the superfluid phonon, and u is the rela- 
tive phase between two counter propagating waves, which 
corresponds to the phonon mode of the stripe. The re- 
quirement of single- valuedness of the wave function cp im- 
poses a constraint that the winding numbers of 9^ and 
0_ must be integers. In terms of 9 and u, that means 
Qv ± Qd are both integers, where gv and qd denote the 
topological charge of vortices (defects in 9) and dislo- 
cations (defects in u), respectively. Three distinct ele- 
mentary topological defects can be identified: (i) a full 
vortex {qv,qd) — (±1,0) (Fig-J^b)) (ii) double disloca- 
tions {qv,qd) — (0,±1) (Fig. [ifc)) and (iii) a half vor- 
tex bound to a single dislocation {q^, qd) — (±1/2, ±1/2) 
(Fig. [l](d)). In two dimensions, proliferation of these 
defects will play crucial roles in finite temperature KT 
transitions. 

For the plane wave phase, superfluid phase is the only 
low-energy degree of freedom, and one can choose the 
ansatz as 
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In this case superfluid vortex is the only type of topolog- 
ical excitation. 



Using these ansatz, one can obtain the effective low- 
energy Hamiltonian following three steps: (i) First 
one can obtain an energy functional J- — {^\Ho — 
g J (PrS'^/2\(p), which characterizes the fluctuation of 9, 
u, ^-1- and ^_ for the stripe phase, or the fluctuation of 
9 and ^ for the plane wave phase, (ii) By expanding it 
to quadratic orders of these fields, we note that both ^+ 
and ^_ (or ^) are massive modes; and (iii) one can then 
integrate out ^+ and ^_ (or ^), which yields an effective 
Hamiltonian for the low-energy physics. The calculations 
are lengthy but quite straightforward, and we present all 
the details in the supplementary material [T^]. Hereafter, 
we shall discuss the results for different cases: 

Stripe Phase with Isotropic SOC: We first consider 
isotropic SOC with rj — 1. In this case, the Hamilto- 
nian has rotational symmetry and the single particle en- 
ergy minimum is in fact a circle in momentum space with 
|k| = K. Nevertheless, the single plane wave state and the 
stripe state are still the only two phases, because the re- 
pulsive density-density interaction {2gQ+g)N^ term does 
not favor density modulation, and therefore rule out a 
superposition state of more than two plane wave states. 
(Unless some other interaction terms are included into 
the Hamiltonian such as for spin-2 bosons [6], which is 
beyond the scope of this work.) 

In this case, the stripe phase can also been called 
"smectic superfluid" in terms of the language of liquid 
crystal physics, because it breaks rotational symmetry 
and translation symmetry along one spatial direction as 
shown in Fig. [ija). Following quite straightforward 
derivations [12j we will arrive at a low-energy effective 
Hamiltonian as 
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where a = ^\ -I- 2k'^ / {mpsg). The most interesting fea- 
ture in "H^J is the absence of (dyu)'^ term, which is a 
manifestation of rotational symmetry. Similar effective 
energy has also been found in the classical smectic liquid 
crystal |13j . A simple argument can be given as follows: 
consider a small fluctuation oi u = qy, i.e. 9-^- = qy and 
9_ — —qy., from the ansatz if, it is clear that this in 
fact corresponds to a change of stripe direction from x to 
Kx+qy, and the module of wave vector is also changed by 
\J (p- ^ K?- — K cx hence, the increase of kinetic energy 
is proportional to q^ instead of . 

In the derivation of the effective theory, it is important 
to keep two massive modes ^_|_ and ^_ first, which rep- 
resent the spin degrees of freedom in the left and right 
moving components, and then integrate them out. Oth- 
erwise, one will end up with a wrong effective Hamil- 
tonian as two decoupled isotropic XY models. This is 
because the rotational symmetry in this system is in fact 
a simultaneous rotation of both spin and space. 

Since the superfluid phase is governed by an 
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FIG. 2: The g — T phase diagram of two-dimensional bosons FIG. 3: The rj — T phase diagram of two-dimensional bosons 
with isotropic SO coupling, g is in the unit of 2k^ /{mps) , with g > 0, where r] = Ky/Kx represents the degree of 
and T is in the unit of Tips /{2m) anisotropy of SO coupling. T is in the unit of Tips/ {2m) 



anisotropic XY model, an algebraic order is expected at 
low temperature. However, the stripe order behaves in a 
very different way due to the absence of the second order 
derivative term along y. The correlation of stripe order 
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Along x-dircction, the integral becomes 
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where q^ = qx\x\ and qy = q-Jy 
direction, the integral becomes 
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where q^ = qx/q^ and qy ^ qy\y 
of stripe "'•^^o'- rio/^Qi™ ov, 



— qx/qy and qy — qy\y\- Hence, the correlation 
order decays exponentially as 
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Eq. ([7| shows that the stripe order becomes disordered 
at any finite temperature. Hence, the broken translation 
symmetry along the stripe direction is restored, and the 
"smectic superfluid" turns into a "nematic superfluid" , 
where only rotational symmetry is broken. Moreover, 
because the relative phase between two components is 
disordered, single boson field can not have broken sym- 
metry, i.e. (-0) = 0. However, boson pair field {i/jtp) 
does not couple to the relative phase, and therefore, still 
exhibits (quasi)-long range order until their superfluid 
phase 29 undergoes KT transition. This phase is named 
"charge-2e nematic superfluid" in the phase diagram Fig. 
[2] where "charge" means particle number. 

This behavior can also been understood from the as- 
pect of topological defects. Because of the absence of 
[dyu)"^ term, the energy of double dislocations (0,±1) 



does not have logarithmic divergence, and fails to com- 
pete with entropy. Therefore dislocations proliferate at 
any finite temperature and melt the stripe order. And 
also, with u field disordered, a (±1/2, ±1/2) defect es- 
sentially becomes a half vortex, and its energy is much 
smaller than a single (±1,0) vortex. Hence, the super- 
fiuid transition of paired bosons is driven by prolifera- 
tion of half vortices. By comparing its energy to entropy, 
one can roughly estimate the transition temperature as 
T'^ = TTps/{8ma), which is reduced by a factor of l/(4a) 
compared to conventional KT transition. 

Plane Wave Phase with Isotropic SOC: Following sim- 
ilar derivations [12 , we arrive at an effective Hamiltonian 
for the plane wave phase as 
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Here the consequence of rotational symmetry breaking is 
now manifested as the absence of (dyd)^ term, which can 
be understood using similar argument presented above. 
Similar analysis can show that the superfluid phase will 
disorder at any finite temperature. Hence, we have com- 
pleted the g — T phase diagram for isotropic SO coupling 
as shown in Fig. [2] 

Stripe Phase with Anisotropic SOC: When 77 < 1, the 
rotational symmetry is broken explicitly. As expected, 
(dyu)'^ no longer vanishes in the low-energy effective 
Hamiltonian given below 
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+idxu)^ + (1 ~ rf){dyuy} , (9) 



and the higher derivative terms can now be ignored. 
However, the resulting XY model is very anisotropic. 
In particular, the coefficient of {dyuY term has to van- 
ish when ry — > 1. For a two-dimension XY model, phase 
transitions are driven by proliferation of different types 
of topological excitations. By calculating the energy of 
different type of topological excitations from the effec- 
tive Hamiltonian Eq. |19[ we can obtain the characteris- 
tic temperatures for proliferation of different topological 
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FIG. 4: The g — rj — T phase diagram of two-dimensional 
bosons with SO couphng. Below the blue surface, on the left 
side of red surface, it is stripe ordered superfluid, and on the 
right side of red surface, it is charge-2e nematic superfluid; 
Below the yellow surface it is plane wave phase, g is in the 
unit of 2k^ / {mps) , and T is in the unit of 7rps/(2m) 



the nature of these two phases, and also will not affect 
the fluctuation analysis discussed above. 

Main predictions of this work can be verified experi- 
mentally. For instance, recently various techniques have 
been applied to study KT transition for spinless bosons; 
Melting of stripe order can be detected through Bragg 
spectroscopy or in-situ image of density. Different topo- 
logical defects can also been differentiated through inter- 
ference technique. Boson pairing may manifest itself in 
the high-order correlation functions. Our results provide 
insights for future experimental study in this new system. 
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It is easy to see that T^j.^ q) is always higher than 
^(±1/2, ±1/2). And for 7] > ?7c = ^/8a'^ /{8a^ + 1), 



T'(±i/2,±i/2) > r'(o,±i)- At ^(o^ii), proliferation of dislo- 
cations first melts the stripe order, but does not destroy 
superfluid phase. Therefore, the system enters "charge- 
2e nematic superfluid" phase. Then, at Ti 



(±l/2,±l/2)> 



half 



vortices proliferate and the system becomes normal. For 
V < Vc, T(±i/2,±i/2) < r(o,±i), (±1/2, ±1/2) defects pro- 
liferate flrst at T'(±i/2,±i/2)5 destroying both superfluid 
and stripe order simultaneously, and driving the system 
into normal phase directly. The rj — T phase diagram in 
the regime g > is shown in Fig. [Sj 

Plane Wave Phase with Anisotropic SOC: For the 
plane wave phase at g < 0, the effective Hamiltonian 
can be derived as 



(11) 



Similarly, the coefficient of (dyO)^ vanishes only when 
77—5-1. In this regime superfluid will turn into normal 
at T = (7rps)-y/l — rf /{2m), and the phase transition is 
driven by integer vortices which is the only type of topo- 
logical excitations in this regime. The transition tem- 
perature is higher than the KT transition driven by half 
vortices in the stripe phase. 

Hence, we have reached a global phase diagram for this 
system as shown in Fig. |4j Our results can be generalized 
to three dimensions. One should also find that stripe 
order melts before superfluid phase is disordered, which 
allows a paired superfluid in three dimensions. Finally, 
we note that, in the real situation, there may be small 
difference between goo and gn, and in some cases there 
is also a Zeeman fleld haz term. They will will not affect 
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Appendix: In this supplementary material, we present some details of derviation of the low-energy effective Hamiltonian. 
For the stripe phase, we have introduced the ansatz 
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By substituting this ansatz into the energy functional and following straightforward calculation, one can obtain 
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(13) 



Here we have dropped some constants and the fast oscillating terms. Expanding this functional to the second order of 9, 
u, = ^+ + f_ and ^2 = ^+ — we obtain (up to some constants) 



J": 



Ps 

8m 



4{V9f+4{Vuf + Uvi.f + ^(V6)' 



+ £ [ + jCl + ri{dye)^2 + Vidyu)^.] + '4^1 (14) 



4^ 4> 

Noting that both and ^2 are massive modes, we can integrate them out and get the energy functional for 9 and u only 
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We can perform expansions in terms of in the denominators. Then we can get 



tie« — tigff + rtefl 



where 



and 



0,9 _ Ps 

2m 



idx9f + (1 



1 + 



lPi^)idy9) 



:^{dy9)V^{dy9) + higher order terms. 



ay" — 



2m 



{dxu}^ + (1 — rf){dyU)^ — ^-^ {dyu)'V^ (dyu) + higher order terms. 
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It is clearly that all four terms {dx9)^ , {dy9)^, (dxuY and {dyuY arc preserved as long as 77 < 1. With all these four 
terms present, other terms involving higher order derivatives of 9 and u (including (dy9)V^{dy6) and (dyu)V^ {dyu)) can 
be dropped in the long wave length limit we are considering. Thus when r? < 1, we end up with an effective Hamiltonian 
for the stripe phase shown as follows: 
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where a = ^Jl + 2k'^ j (rapsg) is defined in the paper. 

In the isotropic limit, i.e. = 1, the terms {dx9)^, (dy9Y still exist, for which all other terms involving higher order 
derivatives of 9 can be ignored. However, {dyu)^ vanishes exactly at the isotropic limit. Therefore, the term (dyu)V^{dyu) 
must be kept. The T-L^g part can be approximated as 
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where an integral by part is performed. Here the second term {dxdyuf can also been ignored in the long wave length 
limit because of the presence (9a; tt)^ term, and one can show that even we include the second term in our calculation, the 
behavior of the correlation functions of u at large distance will not change qualitatively. Considering all of these, we arrive 
at the effective Hamiltonian in the isotropic limit. 
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For the plane wave phase, the ansatz introduced in the paper is 



x+ie j e 



iS/2 



(21) 



(22) 



6 



Similar to the stripe phase, by substituting the ansatz into the energy functional and following straightforward calculations, 
one can obtain 
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Here we have dropped all the constants and fast oscillating terms, and expanded the functional to the second order of 6 
and ^. Again, we find that ^ is a massive field. By integrating out ^, one gets a effective energy functional for 6 only. 
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By the similar argument presented for the stripe phase, only the first two terms in the bracket should be kept, as long as 
Tj < 1. Therefore, our effective Hamiltonian for the plane phase in the anisotropic regime is 



(25) 



In the isotropic limit, the {dy9) vanishes in Eq. 24 The third term in the second line of Eq. 24 becomes important. 
Following the same reason given in discussing stripe phase, we obtain the effective Hamiltonian for the plane wave phase 
in the isotropic limit. 
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